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Abstract 

We study the dynamics of topological defects in the context of "topological inflation" 
proposed by Vilenkin and Linde independently. Analysing the time evolution of planar do- 
main walls and of global monopoles, we find that the defects undergo inflationary expansion 
if i] ~ 0.33mp;, where t] is the vacuum expectation value of the Higgs field and mpi is the 
Planck mass. This result confirms the estimates by Vilenkin and Linde. The critical value 
of r) is independent of the coupling constant A and the initial size of the defect. Even for 
defects with an initial size much greater than the horizon scale, inflation does not occur at 
all if r) is smaller than the critical value. We also examine the effect of gauge fields for static 
monopole solutions and find that the spacetime with a gauge monopole has an attractive 
nature, contrary to the spacetime with a global monopole. It suggests that gauge fields 
affect the onset of inflation. 
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I. INTRODUCTION 



Vilenkin [|T| and Linde [Q] independently pointed out that topological defects can be seeds 
for inflation. The basic idea of this "topological inflation" is as follows. Suppose that we 
have the Higgs field $ a (a = 1, • • • , N), whose potential is 
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V($) = -A($ 2 -v 2 ) 2 , $ = V$ a <5> a , (1.1) 

where rj is the vacuum expectation value and A is a coupling constant. This model gives 
rise to different types of topological defects, i.e., domain walls for N — 1, cosmic strings for 
N = 2 and monopoles for N = 3. Because the center of a defect is in the false vacuum state 
($ = 0), we can expect that inflation occurs near $ = under some condition. Vilenkin 
and Linde claimed that topological defects expand exponentially if and only if 

rj > 0(m Pl ) (1.2) 

for the following reasons. 

(i) If the size of the false vacuum region is greater than the horizon size, it is natural to 
assume that this region undergoes inflationary expansion. The thickness of a domain 
wall in a flat spacetime is given by 

* = A (L3) 



and the horizon size corresponding to the vacuum energy V^(0) is 



Then the condition 5q > H 1 implies rj ~ mpi |i 



* no) 
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(ii) New inflation or chaotic inflation occurs in Friedmann- Robert son- Walker spacetime if 
the slow-rolling condition (|$| 3H\<&\ and $ 2 /2 V(§)) is satisfied. In the model 
( |1.1|) this condition leads to 77 ^> mpij\f^n [0,0. 

(iii) Gravitational effect becomes important if the Schwarzschild radius R g = 2M/m Pl 
is comparable to 5 , where M m (4ti/3)5qV(0). This condition is reduced to r] 



> 



(iv) In the case of a global string or a global monopole there is a deficit angle or a solid 
deficit angle in the static solution. For each defect, this exceeds 2ir or 4-7T if rj > mpi, 
which indicates no static solution with an asymptotically flat region 
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(v) It was shown in that nonsingular static gauge monopole solutions exist only if 
T] ~ mpi, otherwise, the monopole becomes dynamic and may become a Reissner- 
Nordstrom black hole. Linde speculated that in this case the central region of the 
monopole expands exponentially because the condition, r] ~ mpi, is simultaneously a 
condition of inflation (the slow- rolling condition in (ii)). 

In this paper we study topological inflation in detail by numerical analysis. On this 
subject, we need to clarify the following points. 

1. Although the condition for inflation to occur, (|1.2j ), looks plausible, it is based on 
discussions which did not show how the defects really inflate. It is important to verify 
these arguments by numerical analysis which comprehensively includes the gravita- 
tional effect. 

2. Vilenkin and Linde assumed the initial size of defects to be So, which stems from the 
static solution without gravitational effect. We agree that So implies the typical size 
of the defects even in curved spacetime. In the beginning of the phase transition, 
however, the scalar field takes its phase randomly and global distribution may be 
chaotic. Therefore, it is worth studying the dynamics of the defects with various 
initial sizes. Particularly, our interest is the fate of the topological defect when the 
initial size is greater than the horizon scale but rj is not so large. One may conceive 
that such defects also inflate. 

3. The effect of gauge fields is unclear. Linde claimed that gauge fields do not affect 
inflation of strings or monopoles, because they exponentially decrease during inflation. 
We agree that, once inflation occurs, their effect gets smaller and smaller. However, it 
remains an unsettled question whether or not gauge fields affect the onset of inflation. 

In order to answer the first and the second questions, we investigate the evolution of 
planar domain walls and global monopoles in §2. As to the third question, we discuss static 
monopole solutions in §3. §4 is devoted to conclusions. In this paper we use the units 
c = h = l. 

II. EVOLUTION OF PLANAR DOMAIN WALLS AND GLOBAL MONOPOLES 

In what follows we numerically analyse the time-evolution of domain walls (N = 1) in 
a plane symmetric spacetime and of global monopoles (N = 3) in a spherically symmetric 
spacetime. The Einstein-Higgs system is described by the action, 
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(2.1) 



The field equations derived from (2.1) are presented in Appendix. 

For the case of domain walls, we assume that the metric has a form, 

ds 2 = -dt 2 + A 2 (t, \x\)dx 2 + B 2 (t, \x\)(dy 2 + dz 2 ). 



(2.2) 



As an initial configuration of the scalar field, it seems natural to use the same functional form 
as the static solution in a flat spacetime, &fi a t(x) = rj tanh(x / 5) , where S is a parameter of 
the initial width of the domain wall. In a curved spacetime, however, if we consider a single 
domain wall in a vacuum background spacetime, the spacetime becomes closed [[J. There is 
no far region from the wall. For this reason we adopt the periodic boundary condition and 
accordingly change the functional form for each section slightly into 



$(t = Q,x) 



x _ 5f_8_x\ 3 _|_ 3(8_x 



(0 < f < f 
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(2.3) 



The above polynomial has been determined from the following assumptions: (1) 5th-order 
odd polynomial function of x; (2) The first term agrees with that of r]tanh(x/5); (3) It 
connects the constant function $(x) = rj at the point where $' = d^/dx = $" = 0. For 
later convenience, we define a parameter of the initial width of a wall as c = 8/Sq. 

One may wonder if our results really show generic properties of domain walls, because 
the global structure of a spacetime depends on how we choose boundary conditions. We will 
discuss this point in the final section. 

For the case of global monopoles, we assume a spherically symmetric spacetime: 



ds 2 



-dt 2 + A 2 (t, r)dr 2 + B 2 (t, r)r 2 {d9 2 + sin 2 9dtp 2 ). 



For the scalar field, we adopt the hedgehog ansatz: 

<3> a = $(t, r)f a = $(t,r)(sin6 l cosv9,sin6'sin^, cos#). 
As an initial configuration of the scalar field, we suppose 

$(t = 0,r) = $ //at p 



(2.4) 



(2.5) 



(2.6) 



where $fi a t(r) is a static solution in a flat spacetime, which is obtained only numerically, 
and c is an initial size parameter of a monopole normalized by that of the static solution. 

As to the initial value of $ = 9$ / dt, we suppose $ = in both the cases. Following the 
method shown in the Appendix, we solve the constraint equations for setting initial data 
and the dynamical equations for time evolution. 
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Our numerical results are summarized in Figs. 1-5. In Fig.l we show examples of the 
time-evolution of the scalar field, which correspond to a stable domain wall ((a) r\ = 0.2m P i) 
and an inflating wall ((b) rj = O.Qm P i), respectively. The abscissa, X, is defined as a 
proper distance along the x-axis from the origin. The corresponding results for monopoles 
are presented in Fig.2. In the case of (b) we see that the bottle-neck structure appears, 
but we do not find an apparent horizon of a black-hole; this is not a worm-hole. (As for 
an apparent horizon, see the last part of Appendix.) In Fig.3 we draw trajectories of the 
boundaries of a domain wall and a monopole. Here we define the boundary of the defect as 
X b (t) = X or Rb(t) = Br at the position of $ = rj/2. Fig.3 indicates that a domain wall 
or a monopole expands if r\ > OAmpi, while it remains stable if r\ < 0.3mp/. We survey for 
Q.3mpi < rj < OAmpi and 1CT 4 < A < 10 closely, and find the critical value of rj is around 
0.33mp;, regardless of A. This result supports the estimates of Vilenkin and Linde. 

By varying c, we also investigate how the initial size of a domain wall or a monopole in- 
fluences its dynamics. Fig. 4 indicates that the final behavior of a domain wall or a monopole 
is determined only by r\. In the case of (a) or (b) {r\ = 0.2m p ;), even if the initial size is 
much larger than Hq 1 , the configuration eventually approaches to $j/ at (X) or <&fi at (Br). We 
analyse whether or not the cosmological horizon appears or not by searching for an apparent 
horizon, because Hq 1 is not exactly the cosmological horizon. Fig.5 shows that an appar- 
ent horizon really appears and the surface of the shrinking monopole (77 = 0.2mpi,c = 10) 
crosses it. Thus we see that, even for defects with an initial size much greater than the 
horizon scale, inflation does not occur if r\ is smaller than the critical value. 



III. EFFECT OF GAUGE FIELDS FOR STATIC MONOPOLE SOLUTIONS 



In this section we describe the effect of gauge fields for static monopole solutions. We 
consider the SU(2) Einstein- Yang-Mills-Higgs system, which is described by 



S — J d 4 xy/^g 



m pi 
16n 



n 



4 (F;^--(zv£ a ) 2 -vW 



with 



F% = ~ d v Al - ee abc A b ^A c v , D^ a = d^ a + ee abc A b ^ c , 



abc A b 



(3.1) 



(3.2) 



where A a ^ and F® u are the SU(2) Yang-Mills field potential and its field strength, respectively, 
is the covariant derivative. A static and spherically symmetric spacetime is described as 

dg2 = . 2jnjM(E) )e _ MR)dT2 + . 2m- p >M(Ry dR2 + R2(d(j2 + ^ 

(3.3) 
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Under the 't Hooft-Polyakov ansatz, the Yang-Mills potential is written as 

A$ = u?e?*r b l ~ , A a = 0, (3.4) 

where uf is a triad. 

Breitenlohner et al. systematically surveyed regular monopole solutions in this system 
0. According to Fig. 6 in 0, static solutions cease to exist when rj becomes larger than a 
critical value r] cr , which has a little dependence on A/e 2 : 0.20mpi < i] cr (\/e 2 ) < 0.39mp; 
for oo > A/e 2 > 0. (The relations between the parameters in this paper and those in 
are given by T}/mpi = \f^a and A/e 2 = f3 2 /2.) The exact value of i] cr for A/e 2 — > oo 
was analytically obtained 0: i] cr (\/e 2 — > oo) = rripi/\/&K. We also examine this monopole 
system and confirm their results. 

From the two facts of the disappearance of static solutions for large 77 and of the inflation 
of a global monopole for large rj, we might expect that gauge monopoles also inflate for large 
rj. In order to see the effect of gauge fields, we calculate the acceleration of test particles, 
which is governed by the geodesic equation in the coordinate system (|3.3j ): 



d 2 R e 2a ^ d 



dr 2 2 dR 



e 



-2a(R) 



2mpfM(R) 
R 



(3.5) 



where r is the proper time and we have assumed dR/dr = d9/dr = dip/dr = 0. Although 
we cannot discuss the dynamics of a static spacetime by its metric itself, the acceleration 
of test particles gives us some information about the spacetime structure. For example, in 
de Sitter spacetime (M(R) = m 2 Pl H 2 R 3 /2, H = const., a(R) = 0) and in Schwarzschild 
spacetime (M(R) = M = const., a(R) = 0), Q is reduced to d 2 R/dr 2 = H 2 R > and 
d 2 R/dr 2 = —nipfM/R 2 < 0, respectively. The sign of acceleration indicates whether the 
spacetime is repulsive (+) or attractive (— ). In Fig.6 we present examples of the acceleration 
in spacetimes with a global monopole and with a gauge monopole. The figures show that 
the two spacetimes have different properties: the spacetime with a global monopole has a 
repulsive nature 0, while the spacetime with a gauge monopole has an attractive nature. 
Neither sign of acceleration changes even for large 77. Although we have analysed only static 
spacetimes, our results indicates that there may be another possibility that the monopole 
collapses and becomes a black hole. 



IV. DISCUSSIONS 

We have studied the dynamics of topological defects numerically. We have examined 
three questions on "topological inflation" and found the following results. 
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The first subject was simply verifying the arguments of Vilenkin and Linde. Our results 
support their discussions completely and we obtain more exactly the critical value of rj which 
determines whether defects inflate or not. We found that, if 77 ~ 0.33mp; is satisfied, planar 
domain walls and global monopoles inflate. 

The second was to examine the dynamics of defects of various initial size. Our results 
show that in the case of planar domain walls and global monopoles, only rj determines 
whether or not it expands. Even if the initial size of the defect is much greater than the 
horizon scale, it shrinks and approaches a stable configuration if 7] is less than the critical 
value. Some readers may feel this result is surprising, but we can reasonably interpret it as 
follows. First, we found in our analysis that the boundary of a defect (-Xb(t) or Rb(t)) can 
be a spacelike hypersurface, which was also pointed out by Vilenkin [p]]. Thus there is no 
reason that the "horizon scale" prevents a defect from shrinking. Secondly, if rj <C mpi, the 
slow-rolling condition is not satisfied. This condition is expressed only by rj, regardless of A 
and the horizon scale. We therefore conclude that, although the statement "If the size of a 
defect is greater than the horizon scale, it inflates." sounds true, in a strict sense it is not. 

Let us summarize the main points: first, the condition of inflation for domain and for 
global monopoles are the same; secondly, the condition does not depend on the initial size. 
From these two facts, we may understand that whether inflation occurs or not is determined 
not by the global structure of a spacetime but by the local "slow-rolling" condition; our 
results indicate that slow-rolling condition is a necessary and sufficient condition for topo- 
logical inflation. Therefore, although we have assumed a plane symmetric spacetime and the 
periodic boundary condition for a domain wall system, it seems reasonable to conclude that 
our results show generic properties of domain walls. Similarly, we may extend our results to 
the case of global strings. 

The final subject we have investigated is the effect of gauge fields. Comparing the gauge 
monopole static solution with the global monopole one, we find that these spacetimes have 
different characters: one is attractive and the other is repulsive. Therefore, gauge fields 
work to obstruct inflation and their effect cannot be ignored when we discuss the onset of 
inflation. It may be interesting to study time-dependent gauge monopoles, which would give 
us a definite answer to this issue. 
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APPENDIX A: FIELD EQUATIONS AND NUMERICAL METHOD 

In this Appendix we explain how we solve the field equations for planar domain walls 
and for global monopoles. At the end we also summarize how we search for an apparent 
horizon in a global monopole system. 

The variation of ( [2.1| ) with respect to g^ u and $ a yield the Einstein equations, 

8tt 
1 pi 

T, u = d^ a d u $ a - g^-(d a <$> a ) 2 + V($j\, (Al) 



Gnu — T^-llu r\9u.u'^ 9 Tifj/f 

2 m. 



and the scalar field equations, 



□ $« = . (A2) 



(1) Planar domain wall system 

Using the metric fl2.2|) , we write down the field equations ([Al]) and (|A2|) as 



r^- K 2 ( o K 2 ?K) 25" B' 2 2A'B> 8n ( & $' 2 v 
-G = K 2 (3K 2 -2K) -—---—+ — - _ (_ + + F), (A3) 



A 2 B A 2 B 2 A 3 B m£ V 2 2A 2 

1 Gqi = K 2> + ^(3^2 _ R) = *Lq& (A4) 
2 B 

\{G\ + G\ + Gl - Gg) ee K - K\ 2 - 2Kf = ^($ 2 - V) (AS) 

-K 2 -G = K 2 + 1 ^ 2 --K 2 = -j^ Y + — -V) (A6) 

®-K$- — -(-— + —)— + — = A7 
A 2 ^ A B > A 2 d§ 

where we have introduced the extrinsic curvature tensor of t = constant hypersurface, Kij, 
whose components are given by 
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K\ = -j, Kl (= Kt) = -§, (A8) 

and denoted its trace by K = K\. 

In order to set up initial data, we have to solve the Hamiltonian constraint equation (A3) 



and the momentum constraint equation ( |A4|) . We assume the homogeneous and isotropic 
curvature, 

K 

— = K\ = Kl = const., (A9) 
o 



which makes (A4) trivial, and the conformally flat spatial gauge, A = B. Because we have 



adopted the periodic boundary condition such that the period is AS, and assumed reflection 
symmetry, we must keep the condition B'(t,x = 0) = B'(t,x = 25) = 0. By changing K as 
a shooting parameter, we iteratively integrate ( |A"3| ) until the above condition is satisfied. 

Now we have five dynamical variables: A, B, K, K\ and $. (|A8|) , ( |A"5| ) and (|A7|) 
provide the next time-step of A, B, K and $. For the time-evolution of we use (|A6|) 
only at x — 0, and then integrate ( [A"4] ) in the x-direction to obtain other values of K\. In 
this way we have reduced spatial derivatives appearing in the equations, which may become 
seeds for numerical instability. 

(2) Global monopole system 

In a method similar to that applied for the domain wall system, we solve the field 
equations for the monopole system. Here we present the equations and only give comments 
on some differences. Under the assumption of the metric ( p.4| ) and the hedgehog ansatz 
(p.5|), we write down the field equations ( [AT] ) and ( |A2| ) as 



r o _ K 2 ( K 2 K] 2B " B ' 2 2A ' B ' QB ' 2A ' 1 1 
0= 2 ^ 3 2-2 >~^~a^ + ^W~A^B^ + A^~A^ + B^ 

8tt /$ 2 $' 2 $ 2 x /Ain . 

= ^|(T + 2^ + 5V + v )' (A10) 

l Gol = Kl' + (f + l)(3Kt -K) = -^6$'. (All) 

h G \ +G 2 2 + Gl-G°)=k- Kl 2 - 2K? = ^($ 2 - V) (A12) 

■ $" / A' 2B' 2n$' 2$ dV . . > 

— -(-— + — + -) + -^ + — = A13 

A 2 V A B r> A 2 B 2 r 2 d$ v ' 



In this system the homogeneous curvature initial condition ( |A9| ) is not appropriate be- 
cause the far region is asymptotically flat. We thereby suppose A(t = 0, r) = B(t = 0, r) = 1 



S 



and solve the constraint equations (|A10| ) and ( |A11|) to determine Kit = 0,r) and K\{t 



0,r). This treatment is not usually adopted, but it is suitable for this system because we 
obtain 



K 



-K. 



2 



3 2 ~ \ 3m 2 Pi V 2 



» + V), (A14) 



which approaches zero as r increases. The numerical boundary is fixed at r = 105 (5 = cSq). 

In regard to the time-evolution of fTf , we do not have to solve the equation at the origin 
which corresponds to ( |A6|) , because in this case we have the relation K\{t, r = 0) = K(t, r = 
0)/3 from the regularity condition. 

In both the cases (1) and (2), we use a finite difference method with 1000 or 2000 meshes. 
The Hamiltonian constraint equation (|A3| ) or (|A10|) remains unsolved during the evolution 



and is used for checking the numerical accuracy. Through all the calculations the errors are 
always less than a few percent. 

(3) Apparent horizon 

In our analysis for a global monopole system, we investigate the global structure of the 
spacetime, such as the existence of an event horizon of a black-hole, or the existence of a 
cosmological horizon. In numerical relativity, one practical method to see the spacetime 
structure is to search an apparent horizon. The apparent horizon is defined as the outer- 
most/innermost closed 2-surface where the expansion of a null geodesic congruence vanishes. 
For the metric Q2.4Q , the expansion, Q±, is written as 



9 ± = kl. 9 + A; 3 



±;2 ' ^±;3 



K?± {Br) ' 



ABr 



(A15) 



where fc± = (—1, 0, 0) is an outgoing (+) or ingoing (— ) null vector. 

We observe the signs of 0± at all points in the numerical spacetime, and see if there 
is a 2-surface where + or B_ vanishes. In a black hole system, it is proved that, if an 
apparent horizon exists, an event horizon also exists outside (or coincides with) it. Here, as 
well as we search for a black-hole horizon, we use an apparent horizon as a tool for finding a 
cosmological horizon; a cosmological horizon exists inside (or coincides with) the apparent 
horizon. 
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Figure Captions 

Fig.l: Examples of the evolution of a domain wall. We set rj = 0.2mpi in (a), rj = 0.6mpi 
in (b), and c = 1 and A = 0.1 in both cases. The case (a) expresses a stable wall and the 
case (b) expresses an inflating wall. The abscissa is a proper distance along the x-axis from 
the origin. 

Fig. 2: Examples of the evolution of a global monopole. Parameters we choose are the same 
as in Fig.l. In the case of (b) we find that the bottle-neck structure appears. 

Fig. 3: Dependence of the evolution of a domain wall and a monopole on rj. We plot 
trajectories of the positions of $ = rj/2. We set A = 0.1 and c = 1 in both cases, (a) and 
(b) express the dynamics of a domain wall and a global monopole, respectively. A domain 
wall or a monopole expands if rj ~ 0.33mpi. 

Fig. 4: Dependence of the evolution on the initial size. We set A = 0.1 in all cases, (a): 
domain wall, rj = 0.2m Pl . (b): global monopole, rj = 0.2m Pi . (c): domain wall, rj = 0.6mp;. 
(d): global monopole, rj = 0.6mpi. Only rj determines whether or not inflation occurs, 
regardless of the initial size. 

Fig. 5: Trajectories of the apparent (cosmological) horizon and of the boundary of a shrink- 
ing monopole. We set c = 10, rj — 0.2m P i and A = 0.1. The boundary of the shrinking 
monopole crosses the horizon. 

Fig. 6: Acceleration of test particles in spacetimes (a) with a global monopole and (b) with a 
gauge monopole. We set A = 0.1 in (a), A/e 2 = 0.1 in (b), and i] = 0.2mpi in both cases. We 
find that the spacetime with a global monopole has a repulsive nature, while the spacetime 
with a gauge monopole has an attractive nature. 
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